A thermodynamic treatment of the glass transition 
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In undercooled liquids, the anharmonicity of the interatomic potentials causes a volume increase 
of the inherent structures with increasing energy content. In most glass formers, this increase is 
stronger than the vibrational Griineisen volume expansion and dominates the thermal expansion 
of the liquid phase. For a gaussian distribution of inherent states in energy, the generic case, this 
implies a 1/T -temperature dependence of the additional thermal expansion and the additional 
heat capacity at zero pressure. The corresponding compressibility contribution has the Prigogine- 
Defay ratio one. In experiment, one finds a higher Prigogine-Defay ratio, explainable in terms of 
structural volume changes without any energy change. These should always exist, though their 
influence becomes weak in close-packing systems, at the crossover to soft and granular matter. 

PACS numbers: 64.70.Pf, 77.22. Gm 



When an undercooled liquid freezes into a glass at the 
glass temperature T g , the thermal expansion usually de- 
creases by a factor of two to four— _ — . Obviously, the pos- 
sibility to jump from one possible structure to another 
leads to a strong thermal expansion. This possibility gets 
lost as the system enters the glass phase. 

In the glass phase, the thermal volume expansion a g 
has a textbook explanation^ in terms of the Griineisen 
relation for the vibrations 
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where cy g is the heat capacity of the glass at constant 
volume per unit volume, B is the bulk modulus of the 
glass and the Griineisen parameter T describes the vol- 
ume dependence To oc \/V T of the average vibrational 
frequency ZJ. 

The present paper intends the derivation of a similarly 
simple-minded relation for the additional thermal volume 
expansion Aa = a; — a g of the undercooled liquid. The 
additional expansion is ascribed to the volume increase 
of the inherent states with increasing energy content (the 
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FIG. 1: The distribution of inherent states in energy at zero 
pressure. The inset shows the volume distribution of the 
states at E g around the average value V g . 



inherent states are the structurally stable states which 
one obtains when cooling the momentary configuration 
of the liquid to the temperature zcrc 6 ). If the vibrational 
Griineisen expansion is the same for all inherent states, 
one can separate the two influences. As will be seen, this 
simple concept leads to a physical understanding of the 
large Prigogine-Defay ratio at the glass transition 3 ' 4 . 

In first order, the anharmonicity of the interatomic po- 
tential leads to a linear relation between the structural 
energy E and the inherent structure volume V at con- 
stant pressure. Let E g be the average structural energy 
and V g be the average inherent structure volume at the 
glass temperature T g and zero pressure. Then one can 
linearize 



V = V g + a(E- E g ) 



(2) 



for not too large deviations of E from the value E g . The 
anharmonicity factor a, an inverse pressure, is a measure 
for the anharmonicity of the interatomic potential. Note 
that the compressibility of a single inherent state does 
not enter in this relation; it is a glass property, which 
does not require transitions between different inherent 
states to become measurable. 

To get the partition function Z , one has to integrate 
the density g(E) of the inherent states per atom over the 
configurational energy E (see Fig. 1). With the linear 
relation of eq. ([2]) , one then has 



Z = 



g(E) &cp(-/3E(l + pa))dE, 



(3) 



where f3 — 1/kT and the right factor under the integral, 
the Boltzmann factor, contains the pressure product pa. 
One can calculate the average structural energy E per 
atom and the average squared structural energy E 2 per 
atom at zero pressure 
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and 
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E 2 g(E)exp(-PE)dE. 



(5) 



With N atoms in the volume V, the configurational 
part Ac p of the heat capacity at zero pressure per unit 
volume is given by 



Ac,, = — 



N dE 
Vdf 



N ,— =. — 2 N 

(E 2 -E ) 
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(6) 



and since V = V g +a(NE — E g ), the configurational part 
Aa of the thermal volume expansion is 



N dE 
Aa = a VdT =p ' 



(7) 



Finally, the compressibility contribution from the pos- 
sibility to change the average structural energy is easily 
calculated from the derivative of the average energy per 
atom with respect to pressure at constant temperature 
and the pressure zero 



Anpo = a 



N dE Na 2 



V dp VkT 



(E 2 -E 2 ). 



(8) 



The compressibility has the index PD, because it satis- 
fies the Prigogine-Defay relation for a second-order phase 
transition 7 



Ac p Akp D AH 2 AV 2 



(Aa) 2 T (AHAV) 2 



= 1. 



Here AH and AV are the additional enthalpy and vol- 
ume fluctuations from the structural energy changes, re- 
spectively. For completely correlated enthalpy and vol- 
ume fluctuations (an implicit assumption of eq. ([2])), the 
Prigogine-Defay ratio is one. 

According to numerical simulation results^ the density 
g(E) of the inherent states tends to be a gaussian in the 
configurational energy E, with the maximum Eq of the 
gaussian at a value higher than kT g (see Fig. 1) 



g( E ) =3oexp 



-(E-Ep)' 
2^ 



(10) 



where w 2 is the mean squared deviation of E from Eq. 
In this simple case 
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p Vk B T 2 ' 
and the thermal expansion at zero pressure 

Aa 



Vk B T 2 ' 



(11) 



(12) 



Though one cannot expect every glass former to have 
a single gaussian density of inherent states, it turns out 
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FIG. 2: (a) The temperature dependence of the thermal vol- 
ume expansion 9 in B2O3 compared to the sum (continuous 
line) of a constant Griineisen term (dashed line) and the struc- 
tural expansion of eq. ()12p (b) The temperature dependence 
of the heat capacity of the metallic glass vitralloy-1 at ambi- 
(g\ ent pressure^ compared to the sum (continuous line) of the 
vibrational heat capacity of the glass (dashed line) and the 
inherent structure contribution of eq. pip . 



to be easy to find substances which have the 1/T 2 - 
dependence. Fig. 2 shows two examples 9, 1(? , B2O3 and 
the metallic glass vitralloy-1. There are more such cases^. 

But when one measures the Prigogine-Defay ratio at 
the glass transition 2 in B2O3, one finds a value of 4.7, 
much larger than unity. This is contrary to the expec- 
tation of eq. ©. Obviously, there must be additional 
density fluctuations not taken into account by eq. (|8[). 

The inset of Fig. 1 shows the physical reason for these 
additional density fluctuations: The zero pressure vol- 
ume V g at the glass temperature T g is only the average 
of the volume of the inherent states at the zero pressure 
energy E g ; one has to reckon with a mean square devia- 
tion v 2 of their volume distribution. Therefore, eq. ([2]) 
is only valid for the average volume. The volume at con- 
stant structural energy has additional fluctuations, which 
according to the fluctuation-dissipation theorem lead to 
an additional compressibility Ako — v 2 /yfceT. At zero 
pressure, these additional density fluctuations occur at 
constant energy and do neither contribute to the heat ca- 
pacity nor to the thermal expansion. As a consequence, 



one finds the Prigogine-Defay ratio 
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(13) 



where An is the measured value. II is larger than one if 
Ako is larger than zero. 

Consider two inherent states with structural energies 
E\ and E 2 . At T g and zero pressure, they have the 
weighted energy difference (E 2 — Ei)/kT g . After appli- 
cation of a small pressure p, this increases by the factor 
(1 + pa). To restore the weighted energy difference to 
its original value, one has to raise the temperature by 
AT = paTg . Then one returns to the same situation, in 
particular to the same relaxation time r Q of the flow at 
T g . Thus the pressure dependence of T g is given by 



?J±=aT - AaT v. 
dp 9 Ac p 



(14) 



which is in fact the one of the two Ehrcnfcst relations 
for second order phase transitions 11 found to be valid in 
glass formers^. The other one differs by the factor 1/TI 
from the one for a second order phase transition. 

To get an estimate of A^o, consider the structural re- 
laxation processes which bring the shear modulus G of 
the glass down to zero, the elementary processes of the 
flow. After applying a small shear strain to an initial in- 
herent state, the shear relaxation flow processes end up in 
inherent states of the same average structural energy as 
the initial one. The inherent states at higher energy may 
be important for the dynamics as intermediate states, 
but the shear stress release is describable in terms of a 
sum of final transitions at constant average structural en- 
ergy. In general, these transitions will not only change 
the strain state, but the volume as well. Therefore, one 
has to expect a reduction AB of the bulk modulus B of 
the glass by 
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where 5B/SG is the average coupling ratio of the struc- 
tural relaxation processes. Thus 



Ak = 
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B- AB 
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(16) 



so one can determine the average coupling ratio 6B/5G 



from the measurements at T g via 



SB _ BA«(n-l) B 

~5G ~ l + BAfi(n- 1)G' 



(17) 



The question is: Why do some glass formers like B2O3 
have a large Prigogine-Defay ratio (a strong coupling of 
the structural relaxation processes to an external com- 
pression) , while others have a Prigogine-Defay ratio close 
to unity^, i.e. structural relaxations which do not couple 
to an external compression? 



Part of the answer to this question has been given in 
numerical studie o 12 ' 13 of different interatomic potentials. 
These studies have a relaxation time range of nanosec- 
onds, in the best case microseconds. Therefore they do 
not discriminate between vibrations and structural re- 
laxation, but calculate the total enthalpy- volume correla- 
tion. They find a strong correlation in the Lennard- Jones 
potential, a slightly weaker but still strong correlation in 
the MGCU-potential applicable to metallic glasses, but a 
rather weak correlation for hydrogen bonded substances. 

Taking the heavily studied^ Lennard-Jones example, 
the strong enthalpy-volume correlation (equivalent to a 
Prigogine-Defay ratio close to unity) at zero pressure is 
the same as the one for a steep inverse power law po- 
tential with an applied external pressur e ' 1 " ' 1 -'. The in- 
verse power law is i/ r 18 - 9 (7- interatomic distance). This 
shows that the Lennard-Jones potential is rather close 
to the hard-sphere case. The applied external pressure 
needed to hold the atoms together in the equivalent in- 
verse power potential can be estimated from the linear 
potential term of 2.4 r in Lennard-Jones units at the first 
coordination shel l 12 ' 13 ; it is not small. 

The effect of an external pressure on the structural 
states of the inset of Fig. 1, the ones responsible for the 
compressibility Ako, is to change their energies. Thus 
they are no longer at equal energy and begin to contribute 
to the heat capacity and to the thermal expansion. As 
a consequence, the Prigogine-Defay ratio diminishes. Its 
deviation from its zero pressure value should become no- 
table at the critical pressure 
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because then the structural states of the inset of Fig. 1 
have half the energy- volume coefficient of the other ones. 

The strong enthalpy-volume correlation of the 
Lennard-Jones system extends down to low temperatures 
in the glass and even in the crystal^ 2 -, showing that in this 
case the vibrations reflect the properties of the structural 
relaxation. This is consistent with the finding 5B/SG = 
in a Lennard-Jones glass at zero temperature^^. The 
instantaneous affine shear deformation modulus G^ is 
a factor of two higher than the final G, but B^ = B. 
Since B^ and G m have the central-force Poisson ratio 
v = 1/3, this pushes- the final v up to 0.4. The effect 
is due to a non-affine motion of the atoms which lowers 
the shear energy, but does not couple to the compres- 
sion. The non-affine motion is intimately related to the 
boson peak and to the tunneling states which dominate 
the glass behavior at very low temperature s 17 ' 18 as well 
as to the plastic modes responsible for the shear thinning 
in Non-Newtonian flow-iS. 

One can define a (numerically accessible) vibrational 
coupling ratio in the low-temperature glass 

OB Dnn — B . 

(19) 
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which in the Lennard-Jones case is zero, as well as the one 
defined in eq. (|17[) . There seem to be several examples for 



an equality even if the ratio is nonzero^. This supports 
numerical evidence 2 ^ for an intimate relation between the 
soft modes and the structural rearrangements of the un- 
dercooled liquid, a property which glass forming liquids 
seem to share with colloids and granular matter—. 

There are some indications that one has a low SB/SG 
in metallic glasses as well, though their anharmonic 
thermal expansion 23 c\{T g — 0.035 is a factor of ten 
smaller than the Lennard- Jones onei^. Tunneling state 
measurements 24 show an unusual factor of four weaker 
coupling of the tunneling states to longitudinal than to 
transverse waves, consistent with the complete absence of 
a coupling to the compression. Similarly, one does not see 
the Ioffe-Regel limit in the x-ray Brillouin scattering from 
longitudinal waves in a metallic glass 2 - at the boson peak, 
in contrast to measurements in other glass formers^, but 
in agreement with a soft-sphere simulation , which only 
shows the Ioffe-Regel limit for the transverse waves. The 
findings indicate that one has 5B/SG close to zero in 
the rather harmonic metallic glasses as well as in the 
anharmonic Lennard-Jones case, in agreement with the 
numerical finding 12 for the MGCU potential. 

There are two possible reasons for this weak 
compression-relaxation coupling in close packing sub- 



stances: (i) the attractive part of the potential acts as 
a critical pressure in both Lennard-Jones systems and 
metallic glasses (ii) the elementary structural relaxation 
processes in close packing couple only very weakly to the 
compression. The second possibility is supported by two 
atomic models for the structural relaxation in close pack- 
ing, the interstitia l 28 ' 29 and the gliding triangle^, both 
of which couple only to the shear. 

To conclude, the paper presents a thermodynamic de- 
scription of the undercooled liquid which allows to cal- 
culate the additional thermal expansion, the additional 
heat capacity and the additional compressibility above 
T g from the properties of the inherent states. The de- 
scription provides a physical explanation for the large 
measured Prigogine-Defay ratios at the glass transition. 
According to this explanation, a Prigogine-Defay ratio of 
unity is equivalent to a zero coupling of the structural 
relaxation processes to an external compression. In glass 
forming systems, this case is the exception rather than 
the rule. The zero coupling occurs in close packing at 
zero pressure. Inspite of intense numerical studies, the 
reason is not yet fully clear. 

Helpful discussions with Michael Ohl and Herbert 
Schober are gratefully acknowledged. 
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